Mapping spatial persistent large deviations of nonequilibrium surface growth processes 
onto the temporal persistent large deviations of stochastic random walk processes 



O 

o 

(N 

> 
O 

o 

B 

> 
> 



M. Constantin^' ^ and S. Das Sarma^ 

^ Condensed Matter Theory Center, Department of Physics, 
University of Maryland, College Park, Maryland 20742-4111 
^Materials Research Science and Engineering Center, Department of Physics, 
University of Maryland, College Park, Maryland 20742-4111 

Spatial persistent large deviations probability of surface growth processes governed by the 
Edwards- Wilkinson dynamics, Px{x, s), with — 1 < s < 1 is mapped isomorphically onto the tempo- 
ral persistent large deviations probability Pt{t, s) associated with the stochastic Markovian random 
walk problem. We show using numerical simulations that the infinite family of spatial persistent 
large deviations exponents 9x{s) characterizing the power law decay of Px{x,s) agrees, as predicted 
on theoretical grounds by Majumdar and Bray [Phys. Rev. Lett. 86, 3700 (2001)] with the nu- 
merical measurements of 9t{s), the continuous family of exponents characterizing the long time 
power law behavior of Pt{t,s). We also discuss the simulations of the spatial persistence probabil- 
ity corresponding to a discrete model in the Mullins-Herring universality class, where our discrete 
simulations do not agree well with the theoretical predictions perhaps because of the severe finite- 
size corrections which are known to strongly inhibit the manifestation of the asymptotic continuum 
behavior in discrete models involving large values of the dynamical exponent and the associated 
extremely slow convergence to the asymptotic regime. 
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Non-Markovian Gaussian stochastic processes are very 
widely encountered in a large variety of nonequilibrium 
physical problems jj. Considerable theoretical Q and 
experimental ^ 111 efforts have recently been de- 
voted to understanding the first-passage statistics in such 
nonequilibrium systems. Recent work has revealed 
that the time-dependent history of the non-Markovian 
processes can be described by a nontrivial exponent, 9t, 
called the persistence exponent, which depends on the 
dimensionality of the problem, d, and the precise details 
of the non-Markovian nature of the underlying stochas- 
tic dynamics characterizing the phenomenon. The per- 
sistence exponent 9t describes the asymptotic power law 
decay of the persistence probability [Pt{t) oc t~^*] which 
measures the probability that a stochastic variable has 
not changed its characteristics up to time t. As a 
consequence, Ot provides useful quantitative predictions 
concerning the temporal evolution characteristics of a 
given stochastic system. Once the persistence proba- 
bility behavior is found, one can immediately calculate 
the asymptotic behavior of the first-passage probability, 
F(t), which represents the distribution of the time when 
the stochastic variable under consideration first reaches 
a fixed reference value: F{t) = —dPt{t)/dt. In addition, 
one can also obtain the mean first-passage time which 
provides the representative time scale characterizing the 
stability of the dynamical process. Such a time scale 
might be of interest for the study of the evolution of fiuc- 
tuating interfaces or for undestanding the behavior of a 
collection of stochastic spin variables. To be specific, we 
mention that Pt(t) and 9t are the temporal persistence 
probability and exponent, respectively, since we also dis- 
cuss spatial persistence Px (x) and the corresponding ex- 



ponent Ox- 

Of particular interest in the field of surface growth phe- 
nomena is the role played by the dynamics of interfaces 
which are governed by thermal fluctuations. An illus- 
trative category of such interfaces is described by linear 
Langevin equations of the type 

dh{x, t)/dt = -{-V^Y'^h{x, t) + i{x, t), (1) 

where h(x, t) is the step height fluctuation correspond- 
ing to the lateral step position x, at time t, ^{x,t) is 
a white uncorrelated Gaussian noise, and z is the dy- 
namical exponent. It turns out that fluctuating inter- 
faces are of crucial importance at very small scales (i.e. 
nanoscales) involved in the fabrication of current elec- 
tronic devices. In addition to the traditional way of an- 
alyzing various aspects of growth processes based on the 
dynamical scaling behavior of the interface width and 
temporal and spatial correlation functions |M > it has 
been shown that persistence properties provide an ad- 
ditional tool of investigation for understanding the long 
time evolution of surface growth phenomena, due to the 
ability of the nontrivial persistence exponents to iden- 
tify the underlying universality class of the dynamical 
process ^3 ^^^d the presence of the nonlinearities as- 
sociated with the dynamical evolution [1^ . However, 
much broader and more general information can be ex- 
tracted from the natural generalization of the persistence 
through the probability of persistent large deviations 1^ , 
Ptit, s), where — 1 < s < 1. A closely related concept, the 
sign-time distribution, has been introduced in Ref. [isj . 
Pt{t,s) measures the probability that the average sign, 
Sav{t) = Jq sgn[h{x,to + t') — h{x,to)] dt' , remains 
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always above a particular value s up to time t mea- 
sured from the initial time tg. It turns out that Pt{t, s) 
provides an infinite family of temporal persistence ex- 
ponents, Otis), associated with the power law decay of 
Pt{t,s) observed to exist at large time scales for systems 
belonging to different universality classes 0| relevant for 
surface growth. In the limit s — s- 1, the exponent of per- 
sistent large deviations reaches the value of the nontrivial 
persistence exponent, i.e., 9t{s = 1) = Ot- We note that 
the concept of persistent large deviations naturally gen- 
eralizes the concept of persistence exponent from a single 
discrete exponent 9t characterizing the universality class 
to a more general and deeper concept of a continuous 
function, Ot{s), of the persistent exponents characteriz- 
ing the stochastic dynamics. 

The dynamics of the spatially extended systems with 
fluctuations governed by stochastic differential equations 
can be further elucidated by looking, in addition to the 
statistical tools mentioned above, at the spatial analog 
of the temporal persistence, i.e., the spatial persistence 
probability [lallal; Pxix), and its associated exponents. 
Px (x) represents the probability that the height stochas- 
tic variable, measured at a fixed time t, does not reach its 
initial value h{xo,t) up to a longer distance x measured 
from the initial position xq. Theoretical and numeri- 
cal studies [l^ indicate that the power law decay of Px (x) 
is described by two different exponents, the steady-state 
(%s) and the finite initial conditions {9fic) spatial per- 
sistence exponents depending on the selection rules ap- 
plied to Xq : (i) 9ss is obtained if is sampled from the 
entire set of the steady-state configurational sites, and 
(ii) Opjc is obtained if xq is sampled from a subset of 
the steady-state sites characterized by finite height and 
height derivatives. The aim of this paper is to estab- 
lish numerically the concept of spatial persistent large 
deviations probability, Px{x,s) with — l<s<l, asa 
natural generalization of the spatial persistence proba- 
bility concept. We also show that Px{x,s) measured for 
growth processes in the well-studied Edwards- Wilkinson 
|l7| universality class [described by Eq. with z = 2] 
can be mapped isomorphically onto Pt{t, s) of the simple 
random walk stochastic problem. This mapping possibil- 
ity is inspired by the work of Majumdar and Bray (l^, 
who have shown in a recent Letter that the spatial per- 
sistence probability characteristics of growth processes 
involving the interfacial height stochastic variable h{x, t) 
with the dynamics described by Eq. can be mapped 
onto the temporal persistence characteristics of the "ran- 
dom walk" processes of the type d"a;/di" = rj{t), where 
n = {z — d + \)/2 and rj{t) is a white noise as well. The 
purpose of the current paper is to show that this exact 
mapping, as expected, works for the generalized (large 
deviations) persistence probability and the corresponding 
continuous family of persistence exponents as well, and 
to numerically calculate Ox{s) for the important class of 
processes controlled by the Edwards- Wilkinson equation. 



We consider the average sign of the interfacial height 
stochastic variable measured at a fixed time t with re- 
spect to the original value corresponding to the initial 
position Xq, 

SaN(.x) = - sgn[h{xo + x' , t) - h{xQ,t)\ dx' . (2) 
X Jo 

The spatial persistent large deviations probability is de- 
fined, in analogy with its temporal correspondent, as the 
probability that the average sign S'av remains persistently 
above a particular value s, with — l<s<l, uptoa 
longer distance x measured from the initial position xq, 

Pxix, s) = Prob { S'av(a:') > s, Vx' <x}. (3) 

We provide numerical evidence showing that Px{x, s) has 
a power law behavior for x < L, where L is the typical 
length scale in the numerical simulations (i.e. system 
size), independent of the choice of the average sign pa- 
rameter s, 

Px{x,s) ocx~'='^'\ (4) 

where the spatial persistent large deviations exponent 
Ox{s) depends continuously on the parameter s that ap- 
pears in the definition of the probability. The importance 
of Px{x,s) lies in the fact that it provides an infinite 
family of persistence exponents, instead of only one ex- 
ponent as in the case of Px{x). Obviously, Px{x,s — 1) 
and its associated exponent Oxis = 1) are precisely the 
spatial persistence probability and the nontrivial persis- 
tence exponent (^55 or 9pic, depending on the samphng 
procedure applied to xq), respectively. The opposite limit 
s — > —1 is trivial in the sense that Px{x, s = —1) = 1 in- 
dependent of X and as a consequence 0x{s = —1) = 0. 

The s dependence of the temporal persistent large de- 
vations exponents is known exactly for the simple random 
walk case, which is one of the few analytically solved per- 
sistence problems [l^ . 

04 s) = — ^arctanW— ^. 5) 

TT V 1 — S 

The mapping between the temporal properties of 
the random walker (RW) problem and the spatial prop- 
erties of the Edwards- Wilkinson (EW) fluctuating inter- 
faces implies that the expression of Eq. lO also applies to 
the distribution of the spatial persistent large deviations 
exponent as a function of s. This conjecture is verified 
numerically in this study. 

In this paper we have carried out the first application of 
the spatial persistent large deviations concept to the case 
of (l+l)-dimensional fluctuating interfaces characterized 
by the EW dynamical equation. Using the configuration 
of the interface corresponding to a fixed time of the order 
of the time required by the interface width to saturate 
(i.e., t L^), we have computed Px{x^ s) as the fraction 
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FIG. 1: Log-log plot of P{x,s) versus x for the EW equa- 
tion based on the direct numerical integration of Eq. 10 with 
z = 2, using a system of size L = 1000. The average sign 
parameter takes ten different values decreasing from s = 1 
(bottom curve) to s = —0.8 (top curve) with an average sign 
difference As = 0.2. All spatial persistent large deviations 
probabilities show power law decay vs distance for x < L/2. 
The finite-size effects are responsible for the deviations of the 
probabilities from the power law trend at large values of x. 



of lattice sites Xj (with j = 1,2, ...,L — 1) which main- 
tained their stochastic variable 5'av persistently above a 
fixed s value, up to a distance Xj + x. The initial mea- 
surement points, Xj, are sampled from the entire set of 
the steady-state interfacial profile. The numerical inte- 
gration of the stochastic equation is performed using the 
simple forward-time centered-space representation |l9j . 
We have used numerical systems of size L ~ 1000 and we 
have averaged the results over many (~ 1000) indepen- 
dent runs to obtain convergent statistics. 

In Fig. n we show the results for Px{x, s) as a function 
of x for (l-t-l)-dimensional EW interfaces simulated nu- 
merically. We display ten log-log spatial persistent large 
deviations curves versus the distance x for ten values of 
the average sign parameter s (i.e., s = +1, -1-0.8, . . . , — 
0.8). We observe that Rx{x,s) - for x < L/2, 

while for larger values of x and s > there is a down- 
ward deviation of the probability from the power law 
behavior due to finite-size limitations. Except for the 
curve corresponding to s = 1 , which gives the usual spa- 
tial persistence exponent ^55 ~ 0.50, in agreement with 
Refs. mill, all the other curves with s < 1 provide new 
information concerning spatial behavior of the interface 
fluctuations. 

The temporal persistent large deviations probability 
of the random walk model is shown in Fig. [3 We have 
used similar s values, as in the case of Px{x, s) described 
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FIG. 2: Log-log plot of simulated Pt{t,s) vs t for the RW 
problem. The system size is L = 500 and the average sign pa- 
rameter takes ten different values decreasing from s — 1 (bot- 
tom curve) to s = —0.8 (top curve) with As — 0.2 between 
successive probability curves. All temporal persistent large 
deviations probabilities show power law behavior vs time. 

above. Pt{t, s) shows a clear power law behavior versus 
t. We find that Pt{t,s = 1) is characterized by an ex- 
ponent of 0.50, in agreement with the theoretical value 
dt — 1/2. Individual temporal persistent large deviations 
exponents 9t{s) are extracted from the linear regions of 
the log-log plots of Pt{t,s) versus t and they are com- 
pared in Fig. 121 to the corresponding spatial set of expo- 
nents 9rc{s) for the EW interfaces. 

The level of agreement between Px{x,s) correspond- 
ing to the EW dynamical equation and Pt{t,s) of the 
RW case can be seen in Fig.O To generate this figure we 
have used an increment of the average sign parameter (s) 
of 0.1. We observe that the two sets of exponents, Ot{s) 
and Ox{s), overlap very well within the errors of our sim- 
ulations, showing that the mapping procedure involved 
in this study is perfectly applicable. Both cases are in 
agreement with the theoretical prediction of Eq. (|SJ). We 
have also simulated a discrete stochastic growth model, 
the so-called Family model, which is theoretically known 
to exactly belong to the EW universality class [llj. The 
Family model results (not shown here) for Px{x,s) and 
9x{s) are very similar to those shown in Fig. ^ since they 
have identical stochastic dynamics. 

Despite the downward deviation of the probability 
Px{x,s) from the power law behavior due to finite-size 
limitations, we have checked that larger system sizes 
would provide a wider range of distances over which the 
spatial persistence large deviations exponent can be ex- 
tracted with a better precision. This can be seen in Fig.0| 

Another case of interest for epitaxial surface dynam- 
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FIG. 3: dt{s) and 6x{s) vs s as extracted from the power- 
law decay of P{t,s) (for the RW problem) and P{x,s) (for 
the EW fluctuating interfaces), respectively. The increment 
of the average sign parameter is As = 0.1. The continuous 
curve represents the theoretical prediction of Eq. JSJ. 



FIG. 5: Numerical results of Pfic{x) and Pss{x) for the 
LC discrete model with system size L = 200. The mea- 
surements are performed from steady-state configurations. 
Pfic(x) probabilities have been obtained by using three dif- 
ferent band widths, as shown in the legend. Pfic{x) does 
not display a power-law behavior as a function of x over the 
entire range of system size. 
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FIG. 4: Log-log plot of Px{x,a) for s = 1 corresponding to 
the EW equation based on the direct numerical integration 
of Eq. Q with z = 2, using three system sizes, as shown in 
the legend. The straight line represents the fit for L = 10* 
simulation, providing an exponent of 1/2. 



ics is growth under surface diffusion minimizing the local 
curvature, which belongs asymptotically to the MuUins- 
Herring (MH) universality class [i.e. Eq. |^ with 
z = 4] J'he exact mapping prediction by Majumdar and 
Bray 1.51 suggests that the spatial persistence properties 
of the continuum version of the growth models belong- 
ing to this universality class could be mapped onto the 



temporal persistence characteristics of the random ac- 
celeration problem described by the stochastic random 
equation d?x/At^ =r](t) with an analytically known ex- 
ponent of 6t = 1/4 20]. One expects to obtain 0ss = 
and Opic = 1/4 (ISj when measuring the steady-state 
and finite initial conditions regimes of (x), respectively, 
for the MuUins-Herring surface growth dynamics. An ex- 
ample of this case is the (l-l-l)-dimensional model intro- 
duced by Kim and Das Sarma [l^. This discrete solid- 
on-solid atomistic model, the so-called larger curvature 
(LC) model is known to belong asymptotically to the 
MH universality class. As a consequence, we focus on the 
measurement of Ppicix) for the discrete LC model [j^ . 
since Pss{x) is trivially described by a null exponent. 
The definition of Pfic{x) involves the selection of the 
subset of sites characterized by finite height and height 
derivatives. One possibility would be to sample over the 
subset of sites placed on the average level. However, it 
turns out that a system with L = 200, which is the typ- 
ical system size in our simulations, usually has only a 
couple of discrete positions on the average level. For this 
reason, we have sampled over all the lattice sites Xj with 
the height variable (measured with respect to the average 
level) within a band of values characterized by a width 
w [i.e., — w/2 < h{xj) < w/2], where w is taken to be 
smaller than the maximum magnitude of interface fluctu- 
ations. This selection ensured the possibility of sampling 
over a reasonable number of lattice sites presumably suf- 
ficient for good statistical results of Ppic i^) ■ 

In Fig. we show the x dependence of Ppic i^) cor- 
responding to the LC discrete model for three values of 
w: 30, 70 and 110, respectively. The steady-state prob- 
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ability is shown for comparison. We note that Pfic{x) 
does not display the expected power-law behavior as a 
function of x. As the bandwidth w increases, more and 
more lattice sites are included in the sampling subset, 
and Pfic{x) tends to reach the behavior displayed by 
Pss{x). In addition, we observe that when using a nu- 
merical system with L = 200, Pss{x) has a rather linear 
dependence on x, for 50 < x < 200. The impossibility to 
recover the theoretically predicted behavior of Pfic{x) 
may be due to the reduced system size used in our simu- 
lations. This limitation is imposed by the requirement of 
measuring the probability Pfic{x) using an ensemble of 
steady-state configurations that can be achieved only by 
using an extensive computational time ^ L^. We note 
that reducing or increasing the system size by a factor 
of 2 did not produce any qualitative change in the over- 
all behavior of Pfic{x) or Pss{x). In addition, we have 
checked that the direct numerical integration of Eq. 
with z = 4 provides results consistent with the discrete 
LC model. Also, it turns out that similar probability 
curves are obtained for solid-on-solid models belonging 
asymptotically to the molecular beam epitaxy univer- 
sality class (such as the (l-l-l)-dimensional DT model 
[i^V We believe that our problem with the spatial per- 
sistence Px{x) for the LC model belonging to the MH 
universality class arises most likely from the severe 
finite-size problems in simulating systems with large val- 
ues {z = 4) of the dynamical exponent. Large dynami- 
cal exponent implies very slow lateral correlations, which 
considerably complicates studying steady-state behavior 
in the MH universality problem. In fact, this issue is 
very well known in traditional studies of dynamical scal- 
ing involving surface phenomena characterized by a large 
value of the dynamical exponent "24*1 . A large variety of 
stochastic discrete models show long-time transients and 
they cross over very slowly to their corresponding asymp- 
totic behavior. Only extensive simulations of stochastic 
discrete models in the MH universality class can pro- 
vide the asymptotic dynamical scaling associated with 
the continuous limit of the Eq. with z = A. This for- 
bids us from pursuing further measurements of Px{x,s) 
for the MH universality class and checking the validity 
of the mapping procedure, which remains an interesting 
open problem. 

To conclude, we have shown numerically that the spa- 
tial persistent large deviations probability represents a 
possible generalization of the spatial persistence prob- 
ability, providing a useful family of spatial exponents 
for the surface growth phenomena. We have mapped 
these exponents into the family of temporal persistent 
large deviations exponents obtained from the evolution 
of a simple stochastic "random walk" process. We have 
established the validity of this generalization for the 
case of fluctuating interfaces described by the Edwards- 
Wilkinson evolution equation. However, the similar 
problem involving the MuUins-Herring universality class 



remains open since the corresponding discrete LC model 
|22j | simulation shows severe finite-size problems. 
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